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We calculate the quantum corrections to the classical action of a particle with coordinate- 
dependent mass. The result is made self-consistent by a variational approach, thus making it 
applicable to strong-couplings and singular potentials. By including thermal fluctuations, the we 
obtain an effective action whose classical Euler-Lagrange equation describes the motion of a particle 
including quantum and thermal effects. 
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Ch \ I. INTRODUCTION 

_ 

Recently, Jona-Lasinio's group have calculated the first-order quantum correction to the classical equation of 
motion for a particle in one dimension. Both the potential and the kinetic energy terms are modified by terms of 
order h. In this note, we extend these results in four essential ways to make them applicable to nontrivial systems. 
First, we allow the mass to be coordinate-dependent from the outset. This is necessary to describe a larger variety 
I of interesting physical systems, for instance compound nuclei, where the collective Hamiltonian, commonly derived 
■ from a microscopic description via a quantized adiabatic time-dependent Hartree-Fock theory (ATDHF) contains 
coordinate-dependent collective mass parameters. Second, we calculate multi-loop corrections to the effective poten- 
tial. Third, we make the result self-consistent to allow for realistic physical application to systems with strong-coupling 
' or singular potentials. Fourth, we include the effect of thermal fluctuations. 

On the technical side of the one-loop calculation, we differ from jfj] by not using the old-fashioned derivative 
expansion method of Iliopolous et al. 0] , since it cannot be generalized to terms of higher order in derivatives than 
two, but a more powerful method of Fraser which can be extended to any required order. 

II. QUANTUM CORRECTIONS TO EFFECTIVE ACTION 

& . 

Consider a particle with coordinate-dependent mass m(x) moving in a one-dimensional potential V(x). Its 
classical Lagrangian reads 



o 



C(x, x) — \m{x)x 2 ~ V(x), (1) 



where dots indicate time derivatives. According to the path integral formulation of quantum mechanics [g|fj, the 
probability amplitude of the particle initially at position x a at time t a to be found at position Xb at a later time tb is 
given by the path integral 



where A\x] is the classical action 



(x a ,t a \x b ,tb) = I Dxexp 1 -A[x] \ , (2) 



dt£{x,x), (3) 



and the path integral (0) runs over all paths with fixed end points at x(t a ) = x a and x(tb) = Xb- 
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In the path integral formulation, the rules of quantum mechanics appear as a natural generalization of the rules of 
classical statistical mechanics. In statistical mechanics, each volume in phase space is occupied with the Boltzmann 
probability. In the path integral formulation of quantum mechanics, each volume element in the path phase space is 
associated with a pure phase factor exp{i_4[x]/?i}. One may thus consider the quantum-mechanical partition function 



ZQM{tb, t a )[X] 



£>xexp ^ -A[X + x] 
Ti 



(4) 



for the fluctuations x around some background orbit X(t). The result defines an effective action in analogy to the 
free energy in quantum statistics: 



A ctt [X] = -i?ibxZ Q u[X]. 



(5) 



The Euler-Lagrange equations of motion extremizing A cS [X] account for all quantum effects in classical orbits. 

Since Zqm cannot in general be calculated exactly, we must resort to approximations. In the so-called semi-classical 
approximation, A cS is expanded around the classical action, and the quantum corrections are expressed as a series 
expansion in powers of h, also referred to as loop expansion. 

Let x c \(t) be the classical path, solving the classical equation of motion 



m(x c \)x c \ + \m'(x c \)x 2 c \ + V'(x c \) = 0. 
Being an extremum of the action (|^), the first functional derivative of A[x] vanishes at x c \(t): 



SAlx] 



Sx(t) 



dC 



d dC 



d 2 dC 



dx c i dt did dt 2 did 



= 0. 



Hence A[x] has a functional Taylor series starting out like 

S 2 A[x] 



A[x] = A[xd] + \ 



5x{t)5x{t') 



x{t)x{t') 



(6) 



(7) 



(8) 



where x(t) — x(t) — Xd{t) and the quantum-mechanical partition function reads, in the semi-classical approximation, 



■£qm(4, t a )[xd] = exp <{ -^l[x c i] 



Pxexp ( — I dt x(t) 

x a —x b I ^ 



6 2 A[x] 



5x(t)Sx(t) 



i(t) 



(9) 



The path integral in (Q) is Gaussian, and can be calculated analytically (for details see Ref. jg]), yielding an effective 
action 



A ctt [xd] = A[xd] + A^xd] + A ml [xd] 



with the one-loop quantum correction 



Ai[xd] = — ^-Trln 



m{xd)to - V"(xd) - im'{xd)idU ~ m'{xd)id ~ ^"(^cOici 



(10) 



(11) 



where we have defined the time derivative operator Co = — idjdt. A ml contains all multi-loop quantum corrections. 
The functional trace Tr in ( pT| ) contains a time integral J. 6 dt as well as a discrete sum over all eigenvalues cu n of the 
operator Co . Since the summation over discrete eigenvalues introduces unnecessary complications in the calculation, 
we shall replace the sum J2 n /( w n) with an integral J (dw/27r) f(uj). This approximation does not affect the results 
as long as the characteristic frequency along the path 



uj(x) = ^V"(xd)/m(x) 



(12) 



is everywhere much larger than 1/Ai, where At is the inverse time interval where the particle moves very little. The 
precise form of this interval will be specified later. 

If the classical path Xd(t) is taken to be a constant, the determination of the functional trace in ( p"l| ) is straight- 
forward, and the quantum correction to the classical Lagrangian has no explicit time-dependence. We are interested 
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in calculating corrections that are explicitly time-dependent, that is, we look for corrections proportional to x c i(t), 
x c i(t), etc. The time derivative expansion |t| of the correction to the effective action has the general form 



Ai[x d ] = / dt [-V(a;ci(t)) + hZ(x c i(t))x 2 cl + Z 2 {x c i(t))x 2 cl + 



(13) 

The idea behind this derivative expansion Q is to set x c \{t) in ( pd| ) and ( |l5| ) equal to xo+x{t), where xo is a constant, 
and expand both (|ll| ) and (|l^) in powers of x and its derivatives. By comparing the result of both expansions, the 
coefficients V(x), Z(x), etc., are extracted. 

Let us calculate the first two coefficients in (13). An expansion around about Xq up to terms of order x 2 yields 



A\ [x c 



dt 



-V(x ) - V'(x )x - \V'{x Q )x 2 + \Z(x )x + 



(14) 



The corresponding expansion of (O) yields 



Ax = -yTrlnlG- 1 ^) + A(x)\ = -yTr HG" 1 ^)] - yTrln[l + G(w)A(5)], 

where the G _1 (tD) is the inverse free propagator 

G- 1 (u)=m(x )u ) 2 -V"(x ), 

and A(x) contains all time dependent terms: 

A(x) = x[m'(xo)cD 2 — V"'(xq)] + ^x 2 [m" (xo)u> 2 — V""(xq)] ~ ixm! (xq)ui 
— ixxm" (xq)(jj — xm'(xo) — xxm"(xo) — \x m"(xo). 



(15) 



(16) 



(17) 



The first term on the right hand side of (15) can be calculated immediately: the trace is converted into a simple 
integral over time and over the eigenvalues u> of Cj, since V"(xq) is time-independent. Using the integral formula 



do; 



, 2tt [muj 2 - V"(x )] q 
for p — > 0, q — + we obtain 



t -(p+i)(_i)- g [i + (-i)p]r( g - e±i)r(£±^) 

4tt [V " (x ) ] «~ ^ m r [g] 



(18) 



(19) 



Replacing xo by x c \(t), this determines the lowest-order term in the derivative expansion (|13[), in agreement with Ref. 
0. The logarithm in the second term in ( |l5| ) is now expanded up to order x 2 : 



- yTr ln[l + G(u)A(x)] = -yTr[G(w)A(x)] + ^Tr[G(w)A(x)G(c2>)A(x)] 
The first term on the right hand side is calculated using formula (p~q), 



(20) 



$Tr[G(£)A(x)] = / ,1/ - 



m'(zo) 

Tl(xo)] 3 /'- 



■VV"(x ) 



V"'(x ) 



^m(x )V"(x ) 



+ - 



m"(x ) /T7J77 - : V""(xq) 
[m{x )]V^ 1 ° j y/ m (x )V"(x Q ) 



x 2 



1 m"(x ) 
2^/m(x )V"(x y 



(21) 



By comparing the linear term in x with that in (14) we identify 
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V'(zo) = -j 



m'(x Q ) f —— Vi"(xo) 
M^o)] 3 ^ 1 °' ^/m(x )V"(x ) 



(22) 



To find the quadratic terms in (|20|), we must first move all operators ui to the left, and all functions of t to the right. 
Then we can perform the traces independently. To do this, we formulate the commutator [/(i), c2>] = if(t) applied as 
a rule 



f{t)ug{t) = [u+id t ]f{t)g{t), 



(23) 



with the convention that time derivative operators dt act only on the first function on its right. By repeatedly applying 
this rule, the second term in ( p0| ) may be expanded in powers of time derivatives of x, and gives, up to 0(x) 2 

Tt[G(uj)A(x)G(u))A(x)} = Ti[G 2 (i2j)A{x)A(x)} + Tr[G 3 (w) {-2imud t + mdf) A(x)A(x)} 

+ Tr[G 4 (tZ>) (-2imwd t + md 2 t ) 2 A(x)A(x)]. (24) 

In each line, the derivative operators dt act only on the first A(x). After carrying out all integrations, we are left with 



jT±[G(fi)A(x)G(u>)A{x)] 



dt< 



3h[m'(x )} 2 [V"(x )]V 2 , h m'(x )V"'(x ) 



[V"'(x )f 



16 [m(x )] 5 / 2 8 [m{xo)fl' i [V"{x )}y^ 16 [m(x )] 1 /2[^"( a;o )]3/2 

h [V'"(x )} 2 [m{xo)} 1/2 5h m'(x )V'"(xo) 11% [m'(x )} 2 



G4 



[V"{x )f/ 2 



32 [m(a; )] 1 /2[\/»( a;o )]3/2 64 [ m ( a . () )]3/2[y//( a . )]i/2 



f 2 



(25) 



Adding the second term in (|Tj) to the coefficient of the x 2 term in the above expression we obtain —V"(x )/2, 

as necessary for the consistency of the expansion ( |25| ) with ([l3]). From the second term in (p5|), proportional to a; , 
combined to the correspondent term in (21), we extract an xo-dependent contribution to the kinetic energy (|l3|): 



k [V"'{x )] 2 [m{x )Y' 2 5% m'(x )V"'(x Q ) 
32 [V"{x Q )Y>/ 2 



iih 



[m'(x )f 



m"(x ) 



Z(x ) 



We now include all local multi-loop diagrams in the well-known way [g|6J . These smears out the anharmonic part of 
the potential V int {x) = V(x) - V"(x)x 2 /2 to 



16 [m(xo)] 1 /2[T/"( a;o )]3/2 32 [m,{x )} 3 / 2 [V" (x )}^ 2 4 [mixo)] 1 / 2 ^" (x )]V 2 ' 

(26) 



V$\x) 



over the harmonic fluctuation width 



dx' 



y/2ira 2 (x) 



cxp 



i\2 



{x - X 1 ) 
2a 2 (x) 



V int (x') 



= al(x) = (x 2 (t)r = 



h 



2Mlo(x) ' 



(27) 



(28) 



We observe that, for slow enough particle movement, the time interval At determining the range of validity of our 
calculations is given by 



At = 



m(x)al(x) 



Replacing xq everywhere by the position of the background path X, the total effective action becomes 



A cS [X] = 



where the effective potential is given by 



dt 



1 m cS (X)X' z -V°?(X) 



eff/ 



(29) 



(30) 



v cS [x] = ±hu{x) + v£\x). 



(31) 
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calculated within a harmonic path integral containing a trial frequency f2 to be determined by minimization of V^t 
(for details, see Ref. o). The effective mass appearing in Eq. (J3C|) is from (Eq) 



mf{X) = m{X) + 



K[w'{X)f huj'(X)m'(X) 5h 1 [m'(X)} 2 h 1 m"{X) 



16 w 3 (X) 4 uj(X) m(X) 16w(X) m 2 (X) 8w(I) m(X) 



(32) 



If we omit the multi-loop corrections, and the X-dependence of the mass, our effective potential ( p0| ) reduces to the 
expression found by [Q . 

Our result is now made self-consistent as in variational perturbation theory in |^,[|. We simply replace u{X) by a 
trial frequency £l(X) and add to the action the term 



A^ off = 



m(X) 



[u, 2 (X)-tf(X)} a 2 Q (X), 



(33) 



which is the expectation value of the difference between the original and the variational potential. For a slowly 
moving particle, the combined action may be extremized in n(X) for a constant background X(t) = 0. This yields 
the self-consistency equation 



n 2 (x) = 



m{X) d a : 



V a 2{X) = 



1 



m(X) 



dx 2 



(34) 



This is inserted in m^f (X) to obtain the fluctuation-corrected mass term. 

Finally we may include also thermal fluctuations. If the particle moves slowly enough to allow for an approximate 
thermal equilibration at each point of the orbit, we may simply replace the quantum mechanical fluctuation width 
(p8|) at X and Q(X) by the quantum statistical [p~ 



m(X)Q 2 (X) 



nix) a(x) 

' coth ■ v 



2fc B T 



2fc B T 



(35) 



where k^T is the thermal energy. 

The resulting effective action gives rise to a new classical equation of motion, which should help improving our 
understanding of the relation between classical and quantum physics. 
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